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NOMENCLATURE 


Cp  *  pressure  coefficient 

Cpo  *  pressure  coefficient  when  aligned  with  flow 
Cp  «  specific  heat  at  constant  pressure 

m  >  mass  flow  rate 

M  =  Mach  number 

N  =  rotor  speed 

P  “  pressure 

S  *  entropy 

T  »  temperature 

U  “  circumferential  velocity  component 

V  »  absolute  velocity 

W  =  relative  velocity 

X  “  dimensionless  velocity  obtained  by 
dividing  velocity  by  “  s/^CpT^ 

a  *  yaw  angle 

Y  *  ratio  of  specific  heats 

n  -  efficiency 

t<  =  pitch  angle 

j  “  loss  coefficient 


SUBSCRIPTS 

-  Type  A  probe 
»  Type  B  probe 
*  equivalent 

=•  indicator  for  probe  type 
d  “  ideal 

“  relative  frame 
ref  »  referred  quantity 

t  “  total  (or  stagnation)  value 

u  “  circumferential  component 

V  -  absolute  component 

o  =  design  value 

1  *  upstream  of  rotor 

2  »  downstream  of  rotor 


I.  INTRODUCTION 


The  work  presenced  here  reports  and  analyses  test  data  acquired  from  a 
single  stage  axial  flow  compressor.  The  compressor  was  designed  for  super¬ 
sonic  rotor  speeds  with  supersonic  inlet  relative  flow  conditions  over  80%  of 
the  bladespan.  Provisions  were  made  in  the  instrumentation  for  a  variety  of 
data  to  be  acquired.  While  the  full  range  of  compressor  test  speeds  has  not 
been  attempted  as  yet,  the  compressor  has  proved  to  be  a  valuable  tool  for  the 
development  of  new  instrumentation  and  for  the  study  of  transonic  flow  phe¬ 
nomena.  The  testing  so  far  has  been  aimed  at  establishing  the  overall  stage 
performance  map,  measuring  the  rotor  flow  and  developing  a  measurement  tech¬ 
nique  that  resolves  all  three  components  of  the  velocity  vector,  hub-co-tip 
and  blade-to-blade  at  Che  rotor  outlet.  The  latter  is  referred  to  as  Che  Dual 
Probe  Digital  sampling  (DPDS)  technique  and  will  be  described  in  the  course  of 
this  report. 

In  the  course  of  the  test  program,  rotor  through-flow  measurements  showed 
Chat  there  was  a  significant  mismatch  between  the  rotor  and  the  stator  (Ref 
1).  It  was  concluded  that  Che  influence  from  one  on  Che  ocher  would  always  be 
such  that  an  improvement  in  rotor  running  conditions  would  put  the  stator  fur¬ 
ther  off-design  and  vice  versa.  This  was  found  not  to  depend  particularly  on 
speed.  In  order  to  operate  the  rotor  close  to  design,  the  stator  was  removed. 
The  results  given  here  were  derived  from  tests  of  Che  rotor-only  configuration 
at  rotational  speeds  up  to  22,000  RPM. 

The  experimental  program  was  aimed  at  understanding  Che  rotor  aerody¬ 
namics  as  fully  as  possible,  in  particular,  to  determine  the  distribution  and 
magnitude  of  losses.  In  the  present  work,  the  DPDS  technique  was  applied  suc¬ 
cessfully  to  obtain  the  velocity  field  and  the  rotor-relative  loss 
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distribution  everywhere  at  the  rotor  exit 


The  rotor  losses,  hub-to-tip  and 


blade-to-blade ,  were  compared  with  the  predictions  of  loss  models  and  pub¬ 
lished  loss  correlations  for  transonic  compressor  flows. 

II.  EXPERIMENTAL  PROGRAM 

II. 1  TEST  RIG  AND  INSTRUMENTATION 

The  compressor  arrangement  for  the  results  given  here  is  shown  in  Fig  1. 
The  test  rig  is  open  loop.  Ambient  air  enters  the  inlet  settling  chamber 
through  a  microfilter  and  a  throttle  assembly.  It  enters  the  compressor  via  a 
0.45  m  diameter  duct  containing  a  calibrated  flow  nozzle.  From  the  com¬ 
pressor  the  flow  is  turned  and  exhausted  radially  to  atmosphere.  Drive  power 
to  330  KW  is  provided  continuously  by  two  opposed  axial  air  turbine  stages 
powered  by  the  Laboratory  air  supply.  The  rotor  was  designed  for  a  relative 
inlet  Mach  number  of  1.5  at  the  tip,  at  a  speed  of  30,460  RPM.  The  hub-to-tip 
radius  ratio  is  0.5  at  the  inlet  and  0.655  at  the  outlet  with  a  constant 
outer  radius.  The  rotor  design  is  given  in  Ref  2.  The  blading  is 
double-circular-arc  but  with  the  pressure  side  flat. 

The  instrumentation  can  be  divided  into  four  categories;  namely,  fixed 
instrumentation  for  blade  row  performance  measurements,  radial  survey  probes, 
case  wall  static  pressure  tappings  and  high  response  (Kulite)  transducers  with 
associated  computer  software. 

The  blade  row  performance  instrumentation  included  individual  probes  for 
stagnation  temperature  and  pressure  in  the  inlet,  a  distributed  rake  of  20 
kiel  probes  and  3  temperature  probes  in  the  exit  annulus  and  the  flow  nozzle 
measurement  of  flow  rate. 

Four-hole  combination  pressure  and  temperature  probes  (Fig  2)  can  be  tra¬ 
versed  radially  at  measuring  stations  0,  1  and  2.  The  probes  are 


calibrated  as  described  in  Ref  3  to  give  the  time  averaged  velocity  vector,  to¬ 
tal  and  static  pressures.  At  high  inlet  velocities  the  probe  at  station  1  is 
removed  since  the  wake  could  excite  vibrations  in  the  rotor.  The  inlet  condi¬ 
tions  at  station  1  are  then  inferred  from  measurements  at  station  0,  using  the 
now-established  knowledge  of  the  flow  between  these  two  stations. 

Static  pressure  tappings  in  the  compressor  case  wall  cover  a  region  star¬ 
ting  2  chord  lengths  ahead  of  the  leading  edge  and  ending  1.5  chord  lengths 
behind  the  trailing  edge  of  Che  rotor.  Other  Caps  are  located  on  the  compres¬ 
sor  hub  just  behind  the  rotor  trailing  edge. 

Kulite  transducers  used  in  conjunction  with  a  high  speed  digital  data  ac¬ 
quisition  system  provide  two  types  of  measurement,  namely,  unsteady  case  wall 
pressures  and  the  time-resolved  rotor  exit  velocity  vector  distribution.  A 
computer-controlled  sampling  technique  that  allows  the  data  to  be  taken  at 
specific  predetermined  positions  with  respect  to  the  rotor  independent  of 
rotor  speed  is  used  in  both  cases  (Ref  4).  Twelve  wall  pressure  transducers 
are  positioned  between  0.5  chord  length  up  and  down  stream  of  the  rotor.  The 
positions  are  such  chat  there  are  pneumatic  taps  at  the  same  axial  locations, 
but  circumf errentially  displaced.  Blade-to-blade  distributions  of  static 
pressure  are  acquired  with  sufficient  resolution  co  resolve  shock  locations 
within  Che  blade  passages  . 

The  time  resolved  velocity  vector  measurements  involves  two  Kulite  probes 
positioned  in  the  annulus  downstream  of  Che  rotor  and  used  in  a  method  re¬ 
ferred  to  as  the  dual  probe  digital  sampling  (DPDS)  technique. 

II. 2  DPDS  TECHNIQUE 

The  technique  is  described  in  detail  in  Ref  4  and  only  a  brief  summary  of 
the  principles  will  be  given  here.  At  measuring  station  2  (see  Fig  1  and  3) , 


rvo  probes  of  different  design  (Type  A  and  Type  B)  are  mounted  circumferen¬ 
tially  separated  on  the  case  wall  such  that  they  can  be  translated  radially 
and  also  rotated  around  their  tips.  The  probe  outputs  are  sampled  at  the  same 
point  in  the  rotor  frame  of  reference  (at  delayed  times)  with  the  probes  set 
to  nine  angles  with  respect  to  their  axes.  A  procedure  analogous  to  the  cali¬ 
bration  procedure  used  for  multi-sensor  probes  is  applied  to  the  dual  probe 
system  to  derive  liach  number  and  flow  angles  from  the  two  arrays  of  9  measure¬ 
ments  . 

The  Kulite  transducers  in  the  probe  tips  are  protected  by  perforated 
steel  caps  which  enclose  an  internal  volume  of  approximately  0.10  mm^.  In 
contrast  to  the  probes  described  in  Ref  4,  the  present  caps  have  holes  which 
are  oriented  in  straight  lines  across  the  front  face  arranged  parallel  to  the 
probe  shaft  for  the  Type  A  probe  and  perpendicular  to  the  probe  shaft  for  the 
Type  B  probe  (Fig  4).  Due  to  the  different  design  of  the  two  probes,  their 
dependence  on  Mach  number,  pitch  and  yaw  angle  is  quite  different.  In  Fig  5a 
the  output  of  the  Type  A  probe  is  shown  for  one  Mach  number  and  a  range  of 
pitch  and  yaw  angles.  The  probe  is  quite  insensitive  to  pitch  angle,  strongly 
but  symmetrically  dependent  on  yaw  angle.  The  magnitude  of  the  output  in¬ 
creases  with  Mach  number.  Fig  5b  shows  the  Type  3  probe  for  the  same  range  of 
conditions.  The  strong  dependence  on  pitch  angle  shown  was  found  at  all  Mach 
numbers.  Therefore,  while  the  Type  A  probe  is  a  good  indicator  of  the  Mach 
number,  and  independent  of  pitch  angle,  the  Type  B  probe  provides  means  to 
measure  the  pitch  angle.  In  Ref  5  the  use  of  the  probe  data  in  the  calibra¬ 
tion  and  actual  tests  is  described  in  detail.  The  specific  characteristics 
of  these  probes  made  some  changes  necessary.  While  a  Type  3  probe  of  the 
former  design  has  a  clear  maximum  in  its  output  versus  probe  yaw  angle  curve 
when  the  probe  is  aligned  with  the  flow,  this  probe  does  not  have  this 
feature.  Fig  5b  shows  that  for  a  given  Mach  number  the  probe  will  have  a 
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maximum  reading  at  various  yaw  angles  for  different  pitch  angles, 
was  found  that  variations  in  Mach  number  changed  yaw  angle  corresponding  to 
the  maximum  output.  These  variations  did  not  appear  to  be  regular  in  nature. 
In  an  actual  test  the  probe  yaw  angle  corresponding  to  the  maximum  output  of 

the  Type  B  probe  was  taken  as  the  actual  flow  yaw  angle.  This  was  no  longer 

possible  for  the  new  type  probe.  From  the  data  shown  in  Fig  5  only  four  val¬ 
ues  are  required  to  establish  a  calibration  for  the  two  probes  as  system: 

namely,  the  maxima  in  the  output  of  the  Type  A  and  Type  3  probes  and  the  val¬ 

ues  of  the  Type  A  probe  corresponding  to  two  readings  110®  in  yaw  angle  apart. 
These  latter  two  values  are  the  same.  The  center  value  of  the  two  angle  read¬ 
ings  ideally  would  be  the  flow  yaw  angle,  if  the  probe  would  be  symmetrical  in 
yaw  angle.  The  calibration  showed,  chat  the  center  of  Che  110®  yaw  angle 
spread  disagreed  with  Che  flow  yaw  angle  between  1  and  2®.  The  disagreement 
depended  little  on  pitch  angle  but  was  well  behaved  with  Mach  number.  Higher 
speeds  seem  to  improve  the  probe  characteristics.  A  low  order  polynorainal  was 
used  to  express  this  Mach  number  dependence  (Fig  6).  Once  Che  Mach  number  is 
established  in  the  probe's  application,  the  yaw  angle  can  be  corrected  using 
this  polynorainal.  The  four  values  are  treated  like  the  pressure  readings  from 
a  conventional  four  (or  five)  hole  pneumatic  probe.  Two  pressure  difference 
coefficients  are  defined  which  depend  on  Mach  number  and  pitch  angle.  The 
calibration  is  represented  as  surfaces  of  Mach  number  and  pitch  angle  as 
functions  of  Che  two  coefficients,  using  polynomials  in  two  variables  (Ref 
3). 

Another  difference  to  the  former  probes  of  other  tip  geometry  was  ob¬ 
served.  While  the  Type  A  probe  read  the  correct  total  pressure  as  long  as  the 
pitch  angle  did  not  exceed  -5®  and  +15®,  the  new  Type  A  probe  is  more  sensi¬ 
tive  to  pitch  angle  variations.  In  Fig  7  the  dimensionless  pressure  coeffi¬ 
cient  CpoA  averaged  for  all  Mach  numbers  is  plotted  against  flow  pitch  angle. 
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^PoA  is  the  ratio  of  the  difference  between  Type  A  probe  reading  at  an 
aligned  (yaw)  position  and  the  static  pressure  over  the  dynamic  head.  For 
fairly  small  pitch  angles  (8°  or  more)  the  value  of  Cpo^  starts  decreasing 
noticeably.  However,  Cpg^  is  well  behaved  with  pitch  angle  and  is  represented 
by  a  fourth  order  polynominal. 

In  application  in  the  compressor,  the  two  arrays  of  nine  measurements 
(obtained  at  each  point  in  the  rotor  frame  by  rotating  the  two  probes  to  nine 
yaw  angles  varying  from  -60®  to  +80®  about  the  time  averaged  value  read  from 
the  combination  probe),  are  reduced  in  turn  to  obtain  the  local  flow  yaw  an¬ 
gle,  Mach  number  and  pitch  angle.  First,  using  curve  fitting  to  the  Type  A 
probe  data,  the  probe  flow  yaw  angle  is  obtained  as  the  mean  of  the  probe  set¬ 
tings  which  are  calculated  to  give  a  spread  of  110®  between  similar  values  of 
the  probe  output.  The  similar  values  for  the  110®  spread  and  the  maxima  of 
the  Type  A  and  Type  B  probe  outputs  are  also  obtained  by  curve  fitting,  and 
the  two  pressure  difference  coefficients  are  then  calculated.  The  local  Mach 
number  and  pitch  angle  are  obtained  using  the  polynomial  surface  approxima¬ 
tions  to  the  calibration  data.  Using  the  obtained  value  of  Mach  number  the 
flow  yaw  angle  is  corrected  (Fig  6).  The  flow  pitch  angle  is  used  to  deter¬ 
mine  CpoA  which  in  turn  is  used  to  calculate  the  correct  total  pressure  from 
the  Type  A  probe.  The  accuracy  of  the  measurement  is  discussed  in  Appendix  A. 

II. 3  TEST  PROGRAM 

Tests  were  run  at  open  throttle  between  50%  and  70%  of  design  speed, 
blade  row  performance  data  and  wall  Kulite  transducer  data  were  recorded  at 
selected  speeds. 

At  70%  speed,  radial  surveys  of  twenty  measurements  were  made  using  com¬ 
bination  probes  at  stations  0  and  2.  Five  samples  were  recorded  for  each 
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measurement.  At  Che  same  operating  conditions,  the  DPDS  technique  was  applied 
at  five  radial  stations  at  Che  rotor  exit;  namely,  26,  42,  50,  74  and  90  per¬ 
cent  of  Che  blade  span  from  Che  hub.  Data  were  acquired  at  256  positions 
across  two  adjacent  blade  passages.  The  combination  probe  data  were  reduced 
to  give  radial  distributions  of  the  time  averaged  velocity  vector,  stagnation 
pressure  and  temperature.  The  DPDS  data  were  reduced  to  give  blade-to-blade 
distributions  of  velocity  vector  at  the  five  radial  positions.  Distributions 
of  stagnation  pressure  and  temperature,  and  losses  were  also  deduced. 

Additional  tests  were  run  when  periodic  oscillations  at  frequencies  much 
higher  than  blade-passing  frequency  were  observed  on  an  oscilloscope  output 
monitor  of  selected  wall  Kulite  transducers.  The  rotor  speed  was  varied  from 
66%  Co  70%  speed  in  eight  steps  while  the  blade-passing  and  high  oscillation 
frequencies  for  one  transducer  were  obtained  using  a  frequency  analyser  on¬ 
line. 

III.  RESULTS 

The  results  given  here  were  achieved  in  various  test  runs  at  different 
times.  They  are,  however,  presented  in  groups  of  common  subjects  so  that  a 
logical  order  can  be  maintained. 

III.l  TIME-AVERAGED  FLOW 

In  Ref  1  it  is  shown,  that  the  flow  rate  achievable  at  a  full  open  throt¬ 
tle  is  not  sufficient  to  allow  the  correct  incidence  angle  over  Che  whole 
blade  span.  Since  throttling  will  only  reduce  Che  flow  rate,  Che  rotor  in 
flow  would  get  worse.  Thus  no  attempt  has  been  made  to  establish  a  compressor 
performance  map  as  such.  For  a  unthrottled  configuration  the  blade  row  per¬ 
formance  was  measured  at  various  speeds  (Tab  I).  In  order  to  be  able  to  pre¬ 
dict  the  performance  for  higher  speeds  and  to  put  the  measured  performance 
into  perspective  to  individual  speed  lines,  the  performance  map  of  this  rotor 
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was  calculated.  The  through  flow  code  used  was  the  finite  element  Q3DFLO-81 
code  of  Hirsch  (Ref  6  and  7).  Fig  8  shows  the  mesh  used  for  the  calculation 
with  Che  location  of  the  rotor  indicated.  The  performance  map  (Fig  9)  was 
calculated  for  a  speed  range  from  50  to  100%.  Measured  data  (Tab  I)  is  plot¬ 
ted  as  well.  While  the  measured  efficiency  agrees  rather  well  with  Che  calcu¬ 
lated  adiabatic  efficiency,  Che  measured  total  pressure  ratio  is  slightly 
lower  than  the  calculated  one.  For  maximum  efficiency  the  adiabatic  effi¬ 
ciency  does  not  depend  very  much  on  speed,  while  Che  pressure  ratio  changes 
substantially.  Since  Che  speeds  at  which  measurement  data  is  available  do  not 
match  the  calculated  speeds  exactly,  this  is  looked  at  as  a  possible  explana¬ 
tion  for  Che  total  pressure  ratio  discrepancy.  The  measured  flow  rates,  how¬ 
ever,  approximate  those  for  the  calculated  peak  efficiency  rather  well.  Hence 
it  is  assumed  ,  Chat  the  rotor,  if  operated  at  open  throttle,  performs  very 
close  to  its  best  operating  condition.  The  results  of  the  combination  probe 
surveys  were  used  to  derive  the  rotor-relative  flow  conditions,  which  are  re¬ 
quired  for  an  evaluation  of  blade-Co-blade  losses.  Of  particular  importance 
are  the  relative  Mach  number  and  incidence  angle  which  determine  Che  passage 
shock  position  and  strength. 

A  comparison  was  made  with  flow  code  predictions  for  the  given  compressor 
geometry,  test  speed  and  measured  flow  rate.  Results  for  rotor  tn-and  outlet 
conditions  are  shown  in  Fig  10  through  13.  Good  quantitative  agreement  is 
seen  between  measured  and  calculated  inlet  flow  and  good  qualitative  agreement 
is  found  in  the  outlet  flow.  While  the  relative  Mach  number  reached  unity 
only  at  the  tip,  it  exceeded  0.8  over  more  than  50%  of  the  blading. 

The  incidence  angle  to  the  rotor  blading  is  shown  in  Fig  14,  Again,  the 
agreement  between  measurements  and  code  predictions  is  qualitatively  excellent 
and  quantitatively  good.  Also  shown  in  Fig  14  is  the  code  prediction  of  the 
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optimum  (design)  incidence  angle  for  the  given  blading.  Clearly,  in  the  hub 
region  the  incidence  angles  are  larger  than  they  should  be,  for  reasons  dis¬ 
cussed  in  Ref  1  and  8.  Over  the  outer  50%  of  the  blading  the  measured  inci¬ 
dence  angles  were  within  2°  of  design. 

III. 2  DPDS  MEASUREMENTS 

The  DPDS  technique  determines  the  three  dimensional  velocity  vector  in 
the  absolute  frame;  this  includes  a  determination  of  the  pitch  angle,  (which 
is  not  normally  available  from  LDV  or  L2F  measurements).  The  rotor  exit  flow 
from  DPDS  measurements  is  shown  blade-to-blade  and  hub  to  tip  in  Fig  15. 

From  these  quantities  the  rotor  outlet  relative  Mach  number  (Fig  16)  and 
relative  flow  angle  32  (^^8  were  derived.  While  the  Mach  number  varies  as 
is  characteristic  for  a  blade  wake,  the  flow  angle  is  not,  as  assumed  so  far, 
constant  peripherally.  Indeed  rather  large  changes  occur  in  the  blade  wake. 

The  variation  of  pitch  angle  through  the  wake  appears  to  be  quite  large 
and  hence  the  applicability  of  two-dimensional  cascade  flow  concepts  would 
initially  be  suspect.  However,  if  the  pitch  angle  is  calculated  in  the  rela¬ 
tive  system,  it  is  found  to  have  only  about  half  of  the  magnitude  in  the  abso¬ 
lute  system.  This  is  important  if  a  comparison  of  rotor-relative  loss  meas¬ 
urements  with  two-dimensional  steady  cascade  data  is  to  be  meaningful.  The 
rotor-relative  losses  can  be  calculated  from  the  DPDS  measurements  since  the 
DPDS  technique  determines  the  total  pressure  as  well  as  the  velocity  vector. 
The  loss  coefficient  w  for  a  blade  element  is  defined  as  (see  Appendix  3) 


For  the  rotor,  assuming  that  the  inlet  flow  is  axisymmetric,  Pt2Rid>  ^tlR 
?!  are  constant  values  derived  from  combination  probe  data  and  only  Pt2R> 
relative  total  pressure  at  the  rotor  outlet,  varies  circumferentially.  The 
derivation  of  Pt2R  quantities  measured  in  the  absolute  frame  is  given  in 

Appendix  B.  The  calculation  of  Pc2R  carried  out  for  the  DPDS  data 

at  the  five  radii  measured  and  the  result  is  shown  in  Fig  18.  At  26%  span  the 
blade  wakes  are  rather  wide,  probably  due  to  an  incidence  angle  of  6“.  At 
larger  radii  (>55%  span),  the  region  between  two  blade  wakes  ceases  to  be  loss 
free.  At  these  locations,  from  Fig  11  it  can  be  seen  that  the  relative  inlet 
Mach  number  is  0.85  or  larger.  Since  this  exceeds  the  critical  Mach  number, 
it  is  required  that  a  shock  be  present  in  the  passage.  In  Refs  9  and  10  it  is 
proposed  that  such  a  shock  will  produce  a  loss,  that  rises  gradually  from  its 
lowest  value  at  the  pressure  surface  to  its  highest  value  at  the  suction  sur¬ 
face.  In  this  inviscid  model,  in  the  region  of  the  blade  wake,  a  discontin¬ 
uity  between  suction  surface  and  pressure  surface  values  would  occur.  The  in¬ 
crease  in  loss  from  pressure  to  suction  surface  is  largely  the  result  of  the 
increase  in  Mach  number  due  to  the  acceleration  of  the  flow  over  the  suction 
surface.  This  simple  model  appeared  to  explain  the  trends  in  the  measured 
losses  quite  well.  In  order  to  obtain  values  for  the  shock  loss  from  the 
measurements,  the  loss  distribution  outside  the  wake  was  approximated  with  a 
first  order  polynorainal.  The  center  of  the  wake  was  determined  from  the 
blade-to-blade  distribution  of  relative  outlet  Mach  number  (Fig  16)  and  the 
point  where  the  pitch  angle  changed  through  the  time  averaged  value  from  the 
highest  to  the  lowest  value.  Using  a  linear  approximation  the  shock  losses 
were  integrated  at  each  of  the  radii  where  supercritical  conditions  occured. 

Also,  to  obtain  the  total  loss  at  each  radius,  the  relative  outlet  total 
pressure,  Pt2R»  was  mass-averaged  circumferentially  and  used  in  Equ  Cl).  The 
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radial  distributions  of  shock  and  total  losses  derived  from  the  DPDS  tneasure- 


ments,  are  shown  in  Fig  19.  The  total  losses  calculated  from  combination 
probe  surveys  are  also  included  in  Fig  19. 

III. 3  HIGH  RESPONSE  WALL  STATIC  PRESSURE  MEASUREMENTS 


Fig  20.  shows  an  oscilloscope  trace  from  one  high  response  wall  static 
pressure  transducer.  The  transducer  was  located  at  about  mid  chord.  Imposed 
on  the  dominant  blade  passing  frequency  is  an  oscillation  which  varies  in  am¬ 
plitude  depending  on  the  particular  blade  passage.  The  oscillation  occured  at 
a  transducer  and  at  a  blade-to-blade  location  which  suggested  that  it  orig¬ 
inated  from  the  passage  shock.  The  results  obtained  from  the  frequency  ana¬ 
lyzer  showed  that  the  oscillation  frequency  was  some  twenty  times  as  large  as 
the  blade  passing  frequency.  The  variation  of  the  oscillation  frequency  as 
the  rotor  speed  increased  is  shown  in  Fig  21,  It  is  noted  that  as  the  speed 
was  increased,  a  larger  number  of  blade  passages  and  different  transducers 
through  the  rotor  passage  exhibited  the  high  frequency  oscillation.  Also  the 
oscillation  amplitude  appeared  to  increase  with  speed.  While  the  frequency  is 
too  high  for  structual  response  to  be  a  concern,  the  Strouhal  number  based  on 
relative  velocity  magnitude  and  blade  thickness  at  the  tip  is  between  0.1  and 
0.2.  Thus,  shock  oscillation  coupled  with  vortex  shedding  is  a  probable  mech¬ 
anism. 

Assuming  that  the  oscillation  was  caused  by  a  high  frequency  movement  of 
the  passage  shock,  an  approach  proposed  in  Ref  11  and  Ref  12  was  followed  to 
calculate  the  losses  which  would  result  from  the  oscillation.  An  increase  in 
losses  will  result  from  shock  oscillation  since  the  increase  of  entropy  rise 
across  a  shock  is  not  linear  with  increasing  Mach  number  and  thus,  if  the  Mach 
number  varies  about  its  mean  value  (due  to  relative  motion),  the  increase  in 
entropy  rise  for  the  high  relative  Mach  number  will  be  larger  than  the  abso- 


lute  value  of  the  decrease  in  entropy  rise  for  the  low  relative  value.  This 
produces  a  larger  entropy  rise  than  would  be  produced  by  a  steady  shock  at  the 
mean  Mach  number.  Shock  oscillation  amplitude  and  frequency  as  well  as  mean 
value  of  Mach  number  are  therefore  the  parameters  which  determine  the 
magnitude  of  the  losses.  For  the  present  test  conditions,  it  was  found,  that 
due  mainly  to  the  low  inlet  tiach  number  (~1.04)  the  losses  from  the  observed 
shock  oscillation  would  be  negligibly  small. 

IV.  ANALYSIS  OF  THE  LOSS  MEASUREMENTS 

Many  attempts  have  been  made  to  correlate  data  for  losses  measured  in 
compressors  using  models  which  identify  and  sum  various  loss  components.  In¬ 
fluenced  largely  by  the  conclusions  of  Ref  13,  the  losses  obtained  in  the  pre¬ 
sent  study  have  been  compared  with  values  derived  from  the  work  of  Lieblein 
(Ref  14),  Dunker  (Ref  15),  Swan  (Ref  16)  and  Koch  &  Smith  (Ref  17).  It  is  no¬ 
ted  that  the  correlations  used  are  for  cascade  flows.  Since  in  the  compres¬ 
sor,  Che  pitch  angle  relative  to  the  rotor  was  measured  to  vary  peripherally 
no  more  than  +/-7®  through  the  blade  wake,  a  comparison  with  cascade  correla¬ 
tions  is  appropriate. 

The  measurements  resulted  in  data  for  the  shock  loss  and  Che  total  loss 
as  a  function  of  radial  position  (Fig  19).  In  order  to  isolate  shock  loss 
from  Che  total  loss,  an  extrapolation  of  the  inviscid  flow  behavior  Co  a  dis¬ 
continuity  in  Che  blade  wake  was  required.  Clearly,  integration  of  the  vis¬ 
cous  wake  profile  when  shocks  are  present  in  the  outer  flow  gives  a  loss  which 
is  partly  profile  loss  and  partly  shock  loss.  Thus  the  profile  loss  cannot  be 
identified  separately  from  Che  measurements.  The  profile  losses  which  would 
be  expected  from  the  off-design  correlations  of  Lieblein  (Ref  14)  were 
examined.  From  Che  combination  probe  measurements  of  inlet  and  outlet  velo¬ 
cities  and  Che  known  blading  geometry,  the  equivalent  diffusion  factors,  wake 
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niomencum  thicknesses  and  hence  profile  losses  were  calculated  for  each  stream 


surface  following  Ref  14.  Since  the  correlation  for  wake  momentum  thickness 
versus  equivalent  diffusion  factor  involves  an  uncertainty  band,  a  possible 
range  for  the  profile  losses  is  calculated.  The  results  are  shown  in  Fig  22. 
Also  shown  in  Fig  22  are  the  profile  loss  values  which  result  when  two  differ¬ 
ent  through-flow  computer  codes  were  applied  to  the  rotor  geometry.  The 
finite-element  code  Q3DFLO-81  incorporated  the  loss  correlations  of  Koch  & 
Smith  (Ref  17).  The  computer  code  FVVOFF  incorporated  the  profile  loss 
correlations  of  Swan  (Ref  16).  The  latter  results  were  taken  from  Ref  13. 
They  are  seen  to  depart  significantly  from  the  qualitative  variation  expected 
from  the  other  two.  Since  the  center  section  of  the  blade  was  operating 
closest  to  optimum  incidence,  the  higher  profile  losses  given  by  the  FVVOFF 
are  suspect. 

The  shock  losses  deduced  from  the  blade-to-blade  measurements  were  com¬ 
pared  with  the  predictions  of  shock  losses  incorporated  in  the  computer  pro¬ 
grams  Q3DFLO-81  and  FVVOFF.  The  result  is  shown  in  Fig  23.  Q3DFLO-81  uti¬ 
lizes  the  method  of  Koch  &  Smith  (Ref  17)  which  analyses  a  simplified  model  of 
the  flow  incorporating  profile  leading  edge  geometry.  The  Mach  number  and 
leading  edge  thickness  are  the  most  influential  parameters  in  this  model.  The 
model  predicts  no  shock  loss  for  the  present  test  conditions  except  at  the 
blade  tip.  A  more  elaborate  shock  loss  model  developed  by  Dunker  (Ref  15)  is 
incorporated  in  program  FVVOFF.  Similar  to  Moeckel  (Ref  19),  the  model  con¬ 
siders  a  bow  shock  ahead  of  and  a  normal  shock  across  the  blade  passage  to  the 
suction  surface  of  the  next  blade.  The  position  and  strength  of  the  shock, 
and  therefore  the  losses  depend  on  Inlet  Mach  number,  angle  and  rotor  ge¬ 
ometry.  The  agreement  between  the  prediction  with  this  model  and  the  measure¬ 
ments  is  seen  in  Fig  23  to  be  very  good. 
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Since  Che  blade-to-blade  flow  measurements  support  Che  modeling  of  the 
flow  as  being  cascade-lilce  out  to  80%  of  the  span,  it  is  appropriate  to  com¬ 
pare  Che  total  loss  measurements  with  predictions  for  the  sum  of  profile  and 


shock  loss  components.  This  is  shown  in  Fig  24,  where  the  profile  loss  is 
Chat  of  Koch  &  Smith  and  the  shock  loss  Chat  of  Dunker.  The  agreement  is  seen 
Co  be  quite  good.  Whether  similar  agreement  will  be  found  at  higher  transonic 
speeds,  or  whether  stronger  three-dimensional  effects  will  be  found,  will  be 
determined  in  future  tests. 

V.  CONCLUSIONS 

An  axial  compressor  rotor  was  operated  at  moderate  transonic  inlet  condi¬ 
tions  and  the  rotor  flow  field  was  examined  at  inlet  and  outlet.  All  com¬ 
ponents  of  Che  outlet  velocity  vector  were  determined  at  256  points,  blade- 
to-blade,  across  two  identified  blade  passages  at  each  of  5  radial  positions. 
The  losses  in  the  rotor  frame  of  reference  were  derived  as  blade-to-blade  dis¬ 
tributions  and  the  shock  losses  were  evaluated  as  components  of  the  total  loss 
at  each  radius  .  Based  on  the  limited  data  obtained  to  sonic  relative  Mach 
numbers  at  the  tip,  Che  flow  remained  cascade-like  through  the  rotor  without 
strong  three-dimensional  or  case  wall  effects  to  80%  span.  A  combination  of 
the  profile  loss  model  of  Koch  &  Smith  and  shock  loss  model  of  Dunker  was 
found  Co  predict  the  measured  loss  distribution  quite  well.  Measurements  at 
higher  rotor  speeds  are  now  required  Co  examine  the  validity  of  Che  loss 
models  over  a  range  of  Mach  numbers  and  to  investigate  whether  or  not  stronger 
three-dimensional  effects  are  developed.  The  capability  of  the  DPDS  technique 
Co  resolve  Che  radial  as  well  as  tangential  and  axial  velocity  components  of 
Che  velocity.  Including  within  the  rotor  wakes,  was  verified. 
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Figure  4.  View  ot  Kulite  Probe  Tip  Geometry 
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Figure  14.  Flow  Incidence  at  the  Rotor  Leading  Edge 
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Figure  15a.  DPDS  Measurements  of  Blade-to-Blade  Flow 
Mach  Number 
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Figure  15b.  DPDS  Measurements  of  Blade-co-  Blade  Flow 
Yaw  Angle 
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Figure  15c.  DPDS  Measurements  of  31acle-Co-51ade  Flow 
Pitch  .Angle 
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Figure  18.  Measured  Distribution  of  Loss 
Coefficient  31ade-to-Blade 
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Figure  19.  Distribution  of  Losses,  Hub-to-Tip 
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Figure  23.  Shock  Loss  Distribution 
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Figure  26.  Uncertainty  in  Loss  Coefficient  Resulting 
from  Measured  Calibration  Drift 
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Figure  28.  Temperature  Distribution,  Blade-to-Blade 
Derived  from  DPDS  Measurements 


APPENDIX  A 


ACCURACY  OF  THE  RESULTS 
A. 1  COMBINATION  PROBE 

Possible  measuremenc  errors  of  the  combination  probe  could  Influence  Che 
results  In  two  ways.  First,  the  rotor  upstream  conditions,  which  are  used  in 
the  calculation  of  losses,  depend  on  Che  accuracy  of  the  combination  probe. 
Second,  the  XullCe  transducers  used  in  the  DPDS  technique  are  calibrated  on¬ 
line  to  equate  the  time-averaged  measurement  of  Che  dual  probe  system  to  Che 
combination  probe  measurements  of  Che  flow  downstream  of  Che  rotor.  Upstream 
of  Che  rotor  the  combination  probe  la  exposed  to  fairly  low  speeds  In  a  steady 
environment  only.  It  was  shown  in  Che  probe  calibration  process  (Ref  3),  that 
under  these  circumstances  the  probe  will  resolve  Che  velocity  vector  to  better 
than  17,  uncertainty.  However,  in  the  rotor  outlet  plane  the  probe  is  exposed 
to  a  velocity  and  associated  pressure  field,  which  fluctuates  at  blade  passing 
frequency.  In  Ref  20  it  is  demonstrated,  that  a  pneumatic  probe  will  give  an 
Incorrect  time-averaged  pressure  depending  on  the  magnitude  of  pressure  ampli¬ 
tudes  and  frequencies.  For  the  probes  used  here,  it  was  shown  earlier  (Ref 
3),  that  Che  errors  were  negligible  at  rotor  speeds  up  Co  50%  of  design.  The 
evaluation  given  in  Ref  3  was  repeated  for  conditions  found  in  the  present 
work  at  70%  speed  and,  because  of  the  low  amplitudes  in  the  rotor  wakes,  the 
error  was  found  again  to  be  negligible, 

A. 2  VERIFICATION  OF  DPDS  TECHNIQUE  BY  REDUNDANCY 


Where  I  is  either  A  or  B,  depending  on  Che  probe,  the  dependences  of  Cpi  on 
flow  Mach  number,  yaw  and  pitch  angle  are  established  in  the  calibration.  For 
any  given  combination  of  Mach  number  Che  dependence  of  Cpj  on  flow  yaw  angle 
is  a  unique  characteristic  of  the  probe.  In  the  data  reduction  process  only 
four  discrete  values  of  pressure  from  Che  2  probes  are  used.  To  obtain  those 

four  values  however,  data  for  a  range  of  probe  yaw  angle  (9  settings)  is 
acquired,  since  the  local  flow  yaw  angle  is  not  known  prior  to  the 
measurement.  When  Che  data  is  reduced,  values  for  total  and  static  pressure 
are  obtained.  Using  these  and  equation  A(l)  the  pressure  coefficient  Cpj  can 
be  defined  for  each  of  Che  measured  pressures.  Thus  curves  of  Cpi  versus 

probe  yaw  angle  can  be  established  for  each  point  of  measurement  (256  blade- 
to-blade).  If  Che  resulting  curves  match  those  obtained  in  the  steady-flow 
calibration  at  corresponding  Mach  numbers  and  pitch  angles,  Che  measurements 
are  clearly  valid.  Since  only  four  data  values  are  needed  to  derive  the 
velocity  and  flow  angle  magnitudes,  the  demonstration  that  the  remaining  data 
values  from  the  probes  are  consistent  with  the  derived  magnitudes,  provides 
multiply  redundant  verification  of  the  measurements.  In  Fig  25  comparisons 
are  shown  between  measurements  taken  with  the  Type  A  probe  at  90%  span  and 
calibration  curves  at  corresponding  Mach  numbers  and  pitch  angles.  The 
uppermost  figure  shows  results  from  time-averaged  measurements  made  with  the 
dual-probe  system.  (At  each  angle  setting,  the  DC  outputs  of  Che  two  probes 
were  recorded  and  Che  data  reduced  to  obtain  the  time— averaged  velocity  vector 
for  the  complete  rotor,  -  a  measurement  which  is  equivalent  to  chat  made  by 
the  combination  probe).  The  second  figure  is  for  data  taken  at  position 
which  is  in  the  center  of  the  blade  passage.  The  agreement  between  data  and 
the  probe  characteristic  calibration  curve  is  good  in  both 


cases.  The  lower  three  figures  show  results  at  locations  going  into  (’UIZ), 
at  the  center  of  (?Hl6)  and  coming  out  of  (f/138)  the  rotor  wake.  In  contrast 
to  the  results  obtained  with  earlier  probes  (Ref  4),  the  agreement  with  steady 
state  calibration  characteristics  is  quite  acceptable.  The  departures  which 
are  seen  are  isolated  points  and  random  in  character.  The  skewing  of  the 
characteristic  experienced  in  wake  measurements  made  with  the  earlier  probe 
tip  (Ref  4),  was  not  found  with  the  present  geometry  shown  in  Fig  3  and  Fig  4. 

A. 3  CALIBRATION  DRIFT 

In  order  to  account  for  the  temperature  sensitivity  of  the  kulite  trans¬ 
ducers,  Che  calibrations  of  the  Type  A  and  Type  B  probe  transducers  are  per¬ 
formed  on-line  when  test  conditions  have  stabilized.  The  slopes  of  the  trans¬ 
ducers  are  obtained  by  applying  known  pressures  to  the  reference  side.  The 
intercepts  are  derived  by  equating  the  time-averaged  measurements  of  combi- 
,  nation  and  dual  probe  systems.  This  procedure  was  carried  out  before  and 
after  the  blade-to-blade  data  at  nine  angles  were  recorded.  The  average  of 
the  two  calibrations  was  taken  and  the  data  set  was  rejected  if  Che  calibra¬ 
tion  drift  exceeded  2%  in  slope  or  intercept.  The  uncertainty  from  averaging 
the  two  calibrations  is  shown  in  Fig  26  vrtiere  data  are  shown  at  74/i  span  at 
70%  speed.  The  data  were  reduced  using  the  first  and  second  on-line  calibra¬ 
tions.  The  values  of  both  pitch  and  yaw  angle  were  the  same  to  better  chan 
.01®,  the  Mach  number  was  changed  less  than  1%  and  Che  total  pressure  in  Che 
absolute  frame  changed  less  than  1.3%.  While  the  distribution  of  losses 
blade-to-blade  was  little  changed  qualitatively,  the  magnitude  of  the  mass- 
averaged  value  was  increased  or  decreased  about  the  mean  value  by  approxi¬ 
mately  8%.  Similar  changes  were  obtained  in  the  computations  of  Che  shock 
loss.  It  is  the  uncertainty  in  the  loss  values  due  to  the  calibration  un¬ 


certainty  which  is  shown  with  the  data  in  Fig  19 ,  23  and  24 . 


APPENDIX  B 


EVALUATION  OF  PARAMETERS  IN  THE  RELATIVE  FRAME 

In  order  to  compare  the  losses  in  the  rotor  to  those  in  a  two  dimensional 
cascade,  the  losses  in  the  relative  frame  of  the  rotor  must  be  expressed  in 
terms  of  the  measurements  which  are  made  in  the  absolute  frame.  The  losses 
are  defined  as 

5  =  Zlkijd _ 

PtlR  "  Pi 

where  id  denotes  ideal  conditions  as  defined  in  the  T-S  diagram  shown  in  Fig 
27.  The  following  derivations  of  pressures  use  non-dimensional  velocities 
defined  as 

^1-7^  .  Xu2  -  v77  ’  VTJ*  ^2  “ 

where  V{.i  and  Vj.2  are  the  limiting  velocities  given  by 
(Vc  - 

From  Fig  27  we  obtain 

y/y-1  b(3) 

Ptl  Tti 

Tg  -  Tj-i  -  ^  (wf  +  U2  -  Uf  -  vf)  B(4) 

Axial  inflow  results  in  wf  »  Uj  +  V^ 


3(5) 


so  chat 


Pc2Rid  “Pel  (1  + 


B(6) 


PflR  calculated  at  constant  entropy  using 


Pc  IR  _  y/y-1 

Tcl 


•tl 


where 

TclR  -  Tti 

so  that 

is  given  by 
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PclR  -  Ptl  (1  +  4\) 


Pi  -  Pel  (1  -  X^i) 


B(7) 
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In  equations  B(6),  B(7)  and  B(10)  the  measured  (constant)  time  -  average 
quantities  of  Pti,  X^i,  Xy2  are  used.  ^t2id  >  ^tlR  ^1  taken 
to  be  constant  peripherally,  while  Pt2R  varies  blade-to-blade . 

From  Fig  27 
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From  Fig  27 
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and,  rewriting 
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From  the  velocity  triangle 


"■^E  "  Tc2  =•  7^  (^2  -  2V2U2  sin  a2) 


so  that 

=  --V-  (ul  -  2V2U2  sina2)  +  1 
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B(15) 
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Using  Eq  B(ll), 


Pt2R  =  Pt2  (‘^2  -  2  sin  a  Xv2  Xy2  +  D 


Y/'f-l 


B(18) 


To  calculate  ^C2r  from  measurements,  the  individual  readings  of  P-2'  “2 

Xy2  blade-to-blade  have  to  be  taken  from  the  DPDS  data,  while  Xu2  and  Tj^^/T^2 
are  average  quantities. 
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